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1. Introduction 

We study the asymptotic expansions with respect to h of 

P t f(h)±E[A h f(X t )), 

P t f(h)±E[A h f(X t )\& t x ], (1) 
P t f(h)±E[A h f(X t )\X t ], 

with Ahf(Xt) = f(X t +h) — f{X t ), when / : K — > K is a smooth real function, t > is a 
fixed time, X is the solution to the fractional stochastic differential equation 



X t = x+ f b(X s )ds+ f a{X s )dB s , <G[0,T], 
Jo Jo 



(2) 



and JP X is its natural filtration. Here, b, a : M. — > R are real functions belonging to the 
space C£° of all bounded continuous functions having bounded derivatives of all order, 
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while B is a one-dimensional fractional Brownian motion with Hurst index H € (1/2, 1). 
When the integral with respect to B is understood in the Young sense, equation (2) has 
a unique pathwise solution X in the set of processes whose paths are Holder continuous 
of index a E (1 — H,H). Moreover, by, for example, [12], Theorem 4.3, we have, for any 
g:R^Re Cg°, 

g(X t )=g(x) + f g'(X s )a(X s )dB s + f g'(X s )b(X s )ds, t e [0, T]. (3) 
Jo Jo 

The asymptotic expansion of E[/(X^)] with respect to h has been recently studied in 
[1], [7]. In our framework, it turns out to be the case where t = since we obviously have 

E[f(X h )] - f(x) = Pof(h) = Pof(h) = Pof(h). 

In these last references, the authors work in a multidimensional setting and under the 
weaker assumption that the Hurst index H of the fractional Brownian motion B is greater 
than 1/3 (the integral with respect to B is then understood in the rough paths sense of 
Lyons' type for [1] and of Gubinelli's type for [7]). In particular, it is proved in [1], [7] 
that there exists a family 

r = {T 2kH+ i : (k, f)6N 2 , (k, £) ? (0, 0)} 

of differential operators such that for any smooth / : R — > R, wc have the following asymp- 
totic expansion: 

Pof(h) ~ y> 2fcW r 2fe ^(/,cr, &)(*). (4) 

Moreover, in [7], operators T2kH+e are expressed using trees. 

A natural question now arises. Can we also get an expansion of Pt f(h) when t =/= 0? 
Let us first consider the case where B is the standard Brownian motion (which cor- 
responds to the case where H = 1/2). By the Markov property on one hand, we have 
P t f{h) =P t f(h). On the other hand, we always have P t f(h) = E[P t f(h)]. Thus, there 
exist relations between P t f(h), P t f(h) and Ptf(h). Moreover, the asymptotic expansion 
of Ptf(h) can be obtained as a corollary to that of Pof(h) using the conditional expecta- 
tion either with respect to the past of X, or with respect to X t only and the strong 
Markov property 

When H > 1/2, B is not Markovian. The situation regarding P t f(h), P t f(h) and 
Ptf(h) is then completely different and actually more complicated. In particular, we no 
longer have Ptf(h) = Ptf(h) and we cannot deduce the asymptotic expansion of Ptf(h) 
from that of Pof(h). 

The current paper is concerned with the study of possible asymptotic expansions of 
the various quantities Ptf(h), P t f(h) and Ptf(h) when H > 1/2. We will see that some 
nontrivial phenomena appear. More precisely, we will show in Section 3 that P t f(h) 
does not admit an asymptotic expansion in the scale of the fractional powers of h when 
t 7^ 0. Regarding P t f(h), the situations when t = and t > are thus really different. On 
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the other hand, unlike P t f(h), the quantities Ptf(h) and Ptf(h) admit, when t ^ 0, an 
asymptotic expansion in the scale of the fractional powers of h. However, the computation 
of this expansion is more difficult than in the case where t = (as carried out in [1], [7]). 
That is why we prefer to consider only the one-dimensional case. As an illustration, let 
us consider the trivial equation dXt = dBt, t G [0, T], Xq = 0. That is, Xt — Bt for every 
t G [0, T] . We have, thanks to a Taylor expansion, 

fe=l ' k=l 



while by a linear Gaussian regression, when t ^ 
Ptf(h)=E 



l + ^-S + ---)£t + (^ + 



B, 



t 2t 2H 

= HB t f'(B t ) _ Btf'jBt) 2H 
t 2t 2H 

with TV ~ 1) a random variable independent of B t . 

One of the key points of our strategy relies on the use of a Girsanov transformation 
and the Malliavin calculus for fractional Brownian motion. We refer to [3] , [10] for a deep 
insight of this topic. 

We will restrict the exposition of our asymptotic expansions to the case when a = 1. 
Indeed, under the assumption 

(A) the function a is elliptic on R, that is, it satisfies inf^ \o\ > 

and using the change of variable formula (3), equation (2) can be reduced to a diffusion 
Y with a constant diffusion coefficient. 



Y, 



dz 



Moreover, since f Q ^4jy is strictly monotone from R to R under assumption (A), the 
cr-fields generated by X t (resp. by X Sl s <t) and Y t (resp. by Y s , s <t) are the same. 
Consequently, assuming that a = 1 is not at all restrictive since it allows the recovery of 
the general case under assumption (A) . We therefore consider in the sequel that X is the 
unique solution of 

X t =x+ f b(X s )ds + B t , t€[0,T], (5) 
Jo 

with b G Cg° and x G R. 

The paper is organized as follows. In Section 2, we recall some basic facts about 
fractional Brownian motion, the Malliavin calculus and fractional stochastic differential 
equations. In Section 3, we prove that P t f(h) does not admit an asymptotic expansion 
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with respect to the scale of fractional powers of h, up to order n £ N. We eventually 
show, in Section 4, that P t f(h) admits an asymptotic expansion. 

2. Preliminaries 

We first briefly recall some basic facts about stochastic calculus with respect to a frac- 
tional Brownian motion. One may refer to [9], [10] for further details. Let B = {B t )t£[o,T] 
be a fractional Brownian motion with Hurst parameter H £ (1/2,1) defined on a probabil- 
ity space (fi, £/, P). We mean that B is a centered Gaussian process with the covariance 
function E(B s B t ) = iJjj(s,t), where 

RH{s>t)= l {t 2H + s 2H_ lt _ sl 2Hy (6) 

We denote by £ the set of step K- valued functions on [0,T]. Let $) be the Hilbert space 
defined as the closure of £ with respect to the scalar product 

(l[0,i]»l[0,s])j5 =Ri}(t,s). 

We denote by | • |jj the associate norm. The mapping l[o.t] > Bt can be extended to an 
isometry between f) and the Gaussian space Tti(B) associated with B. We denote this 
isomctry tp i— > B(tp). 

The covariance kernel Rn(t,s) introduced in (6) can be written as 

psAt 

R H (t,s)= K H (s,u)K H (t,u)du, 
Jo 

where KH(t,s) is the square-integrable kernel defined, for s < t, by 

K H (t, s) = chs 1 / 2 -" J\u - s) H - 3 / 2 u H -^ 2 du (7) 

with c 2 H = ff(2-2gg-i/2) an< ^ @ ^ ne Beta function. By convention, we set Kfi(t, s) = if 
s>t. 

We define the operator JCh on L 2 ([0,T]) by 

(JC H h)(t)= f K H (t,s)h(s)ds. 
Jo 

Let K.* H : £ — > L 2 ([0, T} ) be the linear operator defined by 

The following equality holds for any <p, ip £ £: 

(<f>, = (K* h 4>,JCh^)v{[o,t\) = E(B(4>)BW). 
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K,* H then provides an isometry between the Hilbert space f) and a closed subspace of 
L 2 ([0,T]). 

The process W = (W t )t£[o,T] defined by 

W t = B((JC* H r 1 (l m )) (8) 
is a Wiener process and the process B has the following integral representation: 

Bt= f K H (t,s)dW s . 
Jo 

Hence, for any <f> £ $), we have 

B(<t>) = wyc* H <i>). 

If b, a £ Cg°, then (2) admits a unique solution X in the set of processes whose paths 
are Holder continuous of index a £ (1 — H, H). Moreover, X has the Doss-Sussman's-type 
representation (see, e.g., [5]) 

X t =<j>(A t ,B t ), t£[0,T], (9) 
with (j) & n d A given, respectively, by 

- — (xi,x 2 ) =cr(</>(xi,X 2 )), (j)(xi,0)=X 1 , Xi,X 2 £M. 
OX 2 



and 



a'(cP(A t ,s))ds)b{0(A t ,B t )), A =x , te[0,T\. 



Let b £ Cg° and X be the solution of (5). Following [11], the fractional version of the 
Girsanov theorem applies and ensures that X is a fractional Brownian motion with Hurst 
parameter H under the new probability Q defined by dQ = rj^ 1 dP, where 

7 ? = cxp^ (^h 1 J o 6WdrJ( S )d^ + ijf (iCh 1 ^b(X r )dr^ (a) da). (10) 

Let be the set of all smooth cylindrical random variables, that is, of the form 
F = f(B(<fii), . . . , B(<p n )), where n > 1, / : R n — > R is a smooth function with compact 
support and <j>i £ f). The Malliavin derivative of F with respect to B is the element of 
h 2 (n,Sj) defined by 

D ° F = Eff" W 1 )' ■ ' s e [0,T]. 
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In particular, D s B t = lm t i(s). As usual, D 1 ' 2 denotes the closure of the set of smooth 
random variables with respect to the norm 

\\F\\l 2 = E[F 2 }+E[\D.F\l]. 

The Malliavin derivative D verifies the following chain rule. If ip : R n — > R is and if 
(i^)*=i,...,n is a sequence of elements of D 1 ' 2 , then ip(Fi , . . . , F n ) £ D 1 ' 2 and we have, for 
any s € [0, T] , 

D s y{F 1 ,...,F n )=Y,-^{Fi,---, F n )D s Fi. 

i=l Xl 

The divergence operator S is the adjoint of the derivative operator D. If a random 
variable u <E L 2 (f2,.f)) belongs to the domain of the divergence operator, that is, if there 
exists c u > such that 

\E(DF,u)^\<c u \\F\\ L 2 for any FeS", 
then S(u) is defined by the duality relationship 

E(F5(u)) = E(DF,u) S} 

for all FeD 1 - 2 . 

3. Study of the asymptotic expansion of P t f(h) 

Recall that P t f(h) is defined by (1), where X is given by (5). 

Definition 1. We say that P t f(h) admits an asymptotic expansion with respect to the 
scale of fractional powers of h, up to order n € N, if there exist some real numbers 
< ol\ < ■ ■ ■ < a n and some random variables Ci,...,C n € L 2 (f2, J^"/ 4 -), not identically- 
zero, such that 

Ptf(h) = dh ai +■■■ + C n h a " + o(h a " ) ash^O, 
where o{h a ) stands for a random variable of the form h a (f>h, with E[0 2 ] — + as h—*0. 

If Ptf(h) admits an asymptotic expansion in the sense of Definition 1, we must, in 
particular, have the existence of a > verifying the following condition: 

lim h~ a P t f(h) exists in L 2 (f2) and is not identically zero. 



However, we have the following. 
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Theorem 1. Let / : K — > R £ Cg° andt£(Q,T]. Assume, moreover, that 

Leb({xeR:/' (x) =0}) = 0. (11) 

Then, as h — > 0, h~ a P t f(h) converges in L 2 (f2) i/ and only if a < H . In this case, the 
limit is zero. 

Remark 1. Since Pof(h) = Pof(h), we refer to Theorem 2 for the case where t = 0. 

Proof of Theorem 1. The proof is divided into two cases. 

(i) First case: a € (0, 1]. Since H > 1/2, let us first remark that h~ a Ptf(h) converges 
in L 2 (ft) if and only if h~ a f '(X t )E[X t + h ~ X t \^ x ] converges in L 2 (f2). Indeed, we use 
a Taylor expansion: 

\f(X t+h ) - f(X t ) - f'(X t )(X t+h -X t )\< Wf'^Xt+H - x t \ 2 , 

so 

\P t f(h) - f(X t )E[X t+h X t \& t x ]\ < |ooE[|X t+h - X t | W]- 
Thus, applying Jensen's formula: 

h- 2a E[\P t f(h) - f'(X t )E[X t+h - X t \^ t x }\ 2 } 

< wni^h-^nnxt+H x t \ 2 \& x ]f 

< l\f"\lh- 2a n\X t+h X t \ A ] = 0(h iH - 2a ). 

Since a < 1 < 2H , we can conclude. 

By (11) and the fact that X t has a positive density on R (see, e.g., [8] Theo- 
rem A), we have that h~ a f'(X t )E[X t +h ~ X t \^ x ] converges in L 2 (f2) if and only if 
h- a E[X t+h - X t \,^ x ] converges in L 2 (f2). 

For X given by (5), we have that ,^ x — . Indeed, one inclusion is obvious, 
while the other can be proven using (9). Moreover, since b £ C£° and a < 1, the term 
h- a m t + b{X s )ds\^ x ] converges in L 2 (D) when h | 0. Therefore, we have, due to 
(5) that h- a V[X t+h ~ X t \^ x ] converges in L 2 (f2) if and only if h- a E[B t+h - B t \^ t B ] 
converges in L 2 (f2). 

Set 

Zf> =h- a E[B t+h -B t \^ t B ] = h- a / (K H {t + h,s)-K H (t,s))dW s , 

Jo 

where the kernel Kjj is given by (7) and the Wiener process W is defined by (8). We 
have 

Var(^ t} ) = h- 2a / (K H (t + h,s)- K H {t, s)) 2 ds 
Jo 
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r-t+h 



r2t>C "l Sl " 2 li (u-s) H - a / 2 u H -^ 2 du^ ds. 



We deduce 

Var(Z^) > h 



C*)N ^ u-2* ( C H \ 2 ,2H-1 



H-l/2, 

f s 1 - 2 "^, + h- sf- 1 ' 2 - (t - s)^ 1 / 2 ) 2 ds 

JO 

-*» f 2 / Yl - ^ ((s + /i)^ 1 / 2 - S ff -V2) 2 ds 



(12) 



H-l/2) J \ t 



with ff (s) = (s + l) ff -V2 _ s ff-i/2. Similarly, 



Var(zW)<fc- !ta ( ] ^ 75 ) a (t + fc) aff - 1 



/* ^-^((t + ft - s)^ 1 / 2 - (t - s)^ 1 / 2 ) 2 ds 
Jo 



-2o 



, x 2H-1 / x 2 

t / \ 1-2H 



/ \ L — zn 



H-1/2J J \~ * 

Note that g 2 (s) ~ (H - |) 2 s 2J? ~ 3 as s — > +oo. So, sg 2 (s) — ► as s-» +oo and 

Jq°° l.9 2 ( s )l ds < +oo since 2H — 3 < —1. Since s i— > sg 2 (s) is bounded on R + , we have, 
by the dominated convergence theorem, that 

1-- -iU 2 (s)ds=/ i L. fl-rWdu 



t J J J u \ h J h 

tends to zero as h — > 0. Thus, 

/•*/?'/ hs\ 1 ~ 2H f°° 
lim / fl-— J g 2 (s)ds = / ,g 2 (s) ds < +oo. 
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Combined with (12)-(13), we now deduce that 

Var(Z<*>) ~ (^^) 2 jT 5 2 (-) as ft - 0. (14) 

If converges in L 2 (fi) as ft — ► 0, then lim^o Var(Z^) exists and is finite. But, thanks 

to (14), we have that lim^o Vai(Z^) = +oo when a > if. Consequently, Zjf does not 
converge in L 2 (f2) as h — » when a > if . 

Conversely, when a < if, we have, from (14), that lim^o Vax(Z^') = 0. Then Z^-^0 
when a < if . 

In order to complete the proof of the first case, it remains to consider the critical case 
where a = if . We first deduce from (14) that Zjf'^+JK(0, a 2 H ), as h — > 0, with 

^ = (^7i) 2 f 92(s)d " 

Let us finally show that the previous limit does not hold in L 2 . Assume for a moment 
that zjf' converges in L 2 (f2) as ft^O. In particular, {Z^ }h>o is Cauchy in L 2 (il). So, 
denoting by Z (t > the limit in L 2 (f2), we have E[Z { e t] Z { s t} ] -> E[|ZW| 2 ] when e,S-> 0. But, 
for any fixed x > 0, we can show by using the same transformations as above that as 
h->0, 

E(^&)-> (h^) 2 ^)=e(I^I 2 ), 

where 

/>OC 

= x g(x u)g(u) du. 
Jo 

Consequently, the function r is constant on ]0, +oo[. The Cauchy-Schwarz inequality 
yields 

\9\h = r(l)=r(V2) = (V2g(2.),g) h2 < \/%(2-)|l= \g\v> = IsIl- 

We thus have an equality in the previous inequality. We deduce that there exists AeM 
such that g(2u) = Xg(u) for all u > 0. Since g(0) = 1, we have A = 1. Consequently, for 
any u > and any integer n, we get 

9(u)=9[^) —> 5(0) = i, 

which is absurd. Therefore, when a~H, zjp does not converge in L 2 (f2) as h — > 0. This 
concludes the proof of the first case. 
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(ii) Second case: a G (l,+oo). If h~ a P t f(h) converges in L 2 (f2), then h~ 1 P t f(h) con- 
verges in L 2 (f2) toward zero. This contradicts the first case, which concludes the proof 
of Theorem 1 . □ 



4. Study of the asymptotic expansion of P t f(h) 

Recall that P t f(h) is defined by (1), where X is given by (5). The main result of this 
section is the first point of the following theorem. 

Theorem 2. Let t G [0,T] and /:l^le Cg°. We write TV for N 2 \ {(0, 0)} . For(p,q)e 
N, set 

J2 P H+q = { (m, n) G N : 2mH + n < 2pH + q} . 



1. If ' t =/= 0, there exists a family {Z^mH+n} (m.n)eN °f random variables measurable 
with respect to X t such that for any (p, q) G Af , 

P t f(h)= Yl Z 2L +n h 2mH+n + o(h 2 ^). (15) 

(m,n)G J 2p H + q 

2. If t = 0, for any (p, q) G Af, we have 

P f(h)=P f(h)=P f(h) 



E f E C ;r / (/,6)( I )]/ l M +"+o(/ i W) 

>.)eJ2pH+ q \/e{04} 2 " l +",l-f|=2m / 



(m,n)£ 

with ci and T/ defined, respectively, by (17) and (19) below. 

Remark 2. 1. In (15), Z^ +1 coincides with the stochastic derivative of X with respect 
to its present t, as defined in [2]. 

2. The expansion (15) allows the expansion of Ptf(h) to obtained for t ^ 0: 

pj(h) = E[p t f(h)} = e n41 H+n }h 2mH+n + o(h 2 ^). 

{m,n)eJ 2 pH+q 

The following subsections are devoted to the proof of Theorem 2. Note that a quicker 
proof of the first assertion seems to be as follows. Once t > is fixed, we write 

rt + h 

X t+h = X t + b(X s )ds + B%\ h>0, (16) 
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where = B t +h — B t is again a fractional Brownian motion. We could then think that 
an expansion for P t f{h) directly follows from the one for Po/(/i), simply by a shift. This 
is unfortunately not the case, due to the fact that the initial value in (16) is not just a 
real number as in the case t = 0, but a random variable. Consequently, the computation 
of E[B { ^\X t ] is not trivial since Bj^ and X t are not independent. 



4.1. Proof of Theorem 2, part (2) 



The proof of this part is actually a direct consequence of Theorem 2.4 in [7]. But, for the 
sake of completeness on one hand and taking into account that we are dealing with the 
onc-dimcnsional case on the other hand, wc give all the details here. Indeed, contrary to 
the multidimensional case, it is easy to compute explicitly the coefficients which appear 
(see Lemma 1 below and compare with [1] Theorem 31 or [7] Proposition 5.4) which also 
has its own interest from our point of view. 

The differential operators T/ appearing in Theorem 2 are recursively defined by 

r (0) (/,&) = 6/', r (1) (/,&) = /' 

and, for I G {0, l} fe , 

r ( i >0 ) (/,&) = &(!>(/,&))', r (Jjl) (/,&) = (r J (/,6)) / (17) 

with (J,0),(/,l)e {0,l} fc+1 . The constants c/ are explained as follows. Set 



dB 



dB u 
dt, 



if i=l, 
if i = 0. 



(18) 



Then for a sequence 



(»!,...,<*)€ {0,1}*, 



wc define 



c 7 = E 



dB 1 



A k [0,l] 



E 



dB\ 



dB 



(ifc-i) 



dB 



(ii) 



(19) 



Set \I\ =X)i<j<fc*i- Equivalently, |/| denotes the number of integrals with respect to 
'dB'. Since B and — B have the same law, note that we have 

C/ = C/ (-l)l 7 l. 

Thus, c/ vanishes when |/| is odd. In general, the computation of the coefficients cj can 
be made as follows. 



*We can also use a rooted trees approach in order to define the Fj's. See [7] for a thorough study, 
even in the multidimensional case and where H > f/3. 
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Lemma 1. Let I € {0, l} k . We denote by J = {ji < • • • < j m } the set of indices j <E 
{1, . . . , k} such that dB^ = dt. We then have that ci is given by 







1 dt Jm - 




10 


Jo 



(fc-im)!(ii-i)! 



n 

k=2 



(jk -Jfc-i)! 



The expectation appearing in the above formula can always be computed using the mo- 
ment generating function of an m-dimensional Gaussian random variable. For instance, 
we have 



E 



E 



E 



1 /-*3 /-*2 



dti dB t2 dB t3 



Uo JO Jo 

1 rt 3 r t 2 



dB tl dt 2 dB t3 



10 Jo Jo 

»1 rts ft 2 



E 



E 



dB tl dB t2 dt 3 

10 JO JO 

pi pt^ pt3 pt 2 

/ / / / dhdhdB^dBt, 
Jo Jo Jo Jo 

/>1 />t4 pta pt-2 

/ / / / dt x dB t ^dHdB ti 

JO Jo JO JO 



dtx dB t2 dB t3 dt 4 




dB tl dt 2 dB ta dt 4 



dB tl dt 2 dt 3 dBt 4 



dB tl dB t2 dt 3 dti 



2{2H + 1)' 

2H-1 
2{2H + l)' 

1 

2(2H + iy 

1 



2(2if +l)(2if + 2) 
H 

(2H+1)(2H + 2Y 
1 

2{2H + 1)(2H + 2) : 
H 

(2H + \){2H + 2Y 

H(2H-1) 
2{2H + 1){2H + 2)'' 

1 

2{2H + l)(2H + 2)' 



(20) 



Lemma 2. When /:l^Ie Cg°, we have 

n—1 p 

f(X h ) = f(x)+J2 E r Ik (fMx) dB Ik {t x ,...,t k ) 



k=U k £{0,l} k 



E / TjMMXt^dB^ih,...^), 

I n £{0,l}" JA ^ ^ 



(21) 



834 



S. Darses and I. Nourdin 



where, again using the convention (18), for g: 



eC 



b ■ 



*A k [0,h] 

Proof. Applying (3) twice, we can write 



f(X h ) = f(x) + f\X s )dB s + (bf')(X s )ds 
Jo Jo 

= fix) + r (1) (/, b)(x) B h + r (0) (/, b) ( x )h 



V / r Ja (/,&)(X tl )dB Ia (t 1 ,t 2 ). 

l{0A}2 JA2[0,h] 



I2G 

Applying (3) repeatedly, we finally obtain (21). 

The remainder can be bounded by the following lemma. 
Lemma 3. If n > 2, e > ( small enough) and g:R^Re C£° are fixed, we have 



□ 



i„e{o,i}» 



A"[0,/i] 



5 (X tl )dS J "(ti,...,i„) 



Proof. This involves a direct application of Theorem 2.2 in [6], combined with the 
Garsia, Rodemich and Rumsey Lemma [4]. □ 

Thus, in order to obtain the asymptotic expansion of P t f(h), Lemmas 2 and 3 say that 
it is sufficient to compute 



E 



/ dB'* fa,... ,t fc ) 

JA k [0M] 



for any Ik € {0,l} fc , with 1 < k < n — 1. By the self-similarity and the stationarity of 
fractional Brownian motion, we have that 



*A k [0,h] 

Hence, it follows that 



/ dB Ik (t u ...,t k )=h H ^ +k -^ f dB Ik {h 

JA k [0,h] JA k [0,l] 



■,t k ). 



E 



dB Ik (h,...,t k ) 



= h H\I k \+k-\I k \ Cik 



I A k [0,h] 

and the proof of point (2) of Theorem 2 is a consequence of Lemmas 2 and 3 above. □ 
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4.2. Proof of Theorem 2, point (1) 
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Let / : R — ► R G C£° and X be the solution of (5). We then know (see Section 2) that X is 
a fractional Brownian motion with Hurst index H under the new probability Q defined 
by dQ = 77- 1 dP, with rj given by (10). Since b : R —> R G Cg°, we observe that 77 G D 1 ' 2 . 
Moreover, the following well-known formula holds for any £ G L 2 (P) n L 2 (Q): 

E[ S |X«] = EQ l*'fl X 'l 



E«h|x,] • 

In particular, 

We now need the following technical lemma. 

Lemma 4. Let £ G B 1,2 (,f)) 6e a random variable. Then for any h > 0, the conditional 
expectation E[C,{f {B t +h) — f (B t ))\B t ] is equal to 

t+h 



pi pi-t-n 

H(2H-l)f'(Bt)J duE[D u C\B t ]J^ \v-u\ 2H - 2 dv 



+ -t- 2H f'(B t )(B t E[C\B t } ~ E[(D(, l m ) s ,\B t ])(h 2H - (t + hf H + t 2H ) 
- ^t- 2H f"(B t )E[C\B t ] 

rt+h 

x J {{v-t) 2H - 1 -v 2H - l ){t 2H +v 2H ~(v-t) 2H )dv 

+ ±f"(B t )E{C\Bt]((t + h) 2H -t 2H ) 
rt+h 



pi pt+n 

+ H(2H-1) du \v-u\ 2H - 2 E[D u ({f(B v )-f(B t ))\B t ]dv 
Jo Jt 

+ Ht- 2H B t j' + {{v ~ t) 2H ~ l - ^ 1 )E[C(/'(fi„) - f\B t ))\B t ] dv (22) 

rt+h 

Ht~~ 2H J ((v t) 2 *- 1 v 2H -i)E[{D(, l [0)t] > fl (/'(B„) - f(B t ))\B t ] dv 



H 



t+h 



_ ^ t ~ 2H I ((v - t) 2 "- 1 - v 2H - l ){t 2H + v 2H -(v- t) 2H ) 
xE[af"(B v )-f"(B t ))\B t ]dv 
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pt+h 

+hj^ Eiccr^)-/"^))!^]^- 1 ^. 

Proof. Let j:R^IeCj. We can write, using the Ito formula ([9], (5.44), page 294) 
and basics identities of Malliavin calculus, 

E[{g(B t )(f(B t+h )-f(B t ))] 

pt+h 

= E[(g(B t )6(f(B.)l [tit+h] )}+H J E[Cg{B t )f" {B^v 2 *- 1 dv 

-EL 9 (i?0(^C,l[M + h]/ , (i?.)) fl ]+E[C5'(^)(l[o,t],l[ t , t+ / l ]/'(B.)) fi ] 
pt+h 

+ Hj EKgiB^f'iB^v^dv 

pT pt+h 

= H(2H-1) du \v-u\ 2H - 2 E[f'{B v )g[B t )D u C]dv 
Jo Jt 

pt pt+h 

+ H(2H-1) du (v-u) 2H - 2 E[{f(B v )g'(B t )]dv 
Jo Jt 

+ H £ + ' EiCgiB^f'iB^v 211 - 1 dv. 



But. 



E[Cg(B t )f'(B v )B t ] = E[{D(Cg(B t )f'(B v )),l m )^} 

= E[g(B t )f'(B v )(D(,l [0tt] )z}+El(g'(B t )f'(B v )]t 2H 
+ \E[Cg{B t )f"{B v )]{t 2H + v 2H -(v- t) 2H ). 

Consequently, 

E[(g(B t )(f(B t+h )-f(B t ))} 

pT pt+h 

= H(2H-1) du \v-u\ 2H - 2 E[f(B v )g(B t )D u (}dv 
Jo Jt 

pt pt+h 

+ H(2H-l)tr 2H du (v-u) 2H - 2 E[(g(B t )f'(B v )B t ]dv 
Jo Jt 

-H(2H-l)r 2H f du [ t+h (v~u) 2H - 2 Elg(B t )f(B v )(D(,l [Q . t] )^}dv 
Jo Jt 

1 /"* pt+h 

- -H(2H - l)t- 2H J du J (v - u) 2H - 2 El(g(B t )f'(B v )} 



x (t 2H +v 2H -{v-t) 2H )dv 
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rt+h 



/t+n 
nCg(Bt)f"{B v )]v 2H -Uv. 



We deduce that 

nCU(B t+h )-f{B t ))\B t ] 

r T r t+h 



H(2H-1)J du£ \v-u\ 2H ~ 2 E[f'{B v )D u C,\B t ]dv 

rt+h 

+ Hr 2H B t J ((« - tf H ~ x - v^-^EK f'(B v )\B t ] dv 
Ht- 2H / <+ ((« - tf H ~ l - v 2H - l )E[f(B v )(DC, l [0> t]> a |B t ] dv 



H 



t+h 



t~ 2H / ((v-t) 2 "- 1 -v 2H - v )(t 2H +v 2H -(v~t) 2H )E{(f"(B v )\B t }dv 
* Jt 

+ H J*^ E[C/ / '(S„)lS t ]« Mr - 1 d». 
Finally, (22) follows. □ 

First, we apply the previous lemma with £ = 77, E = and B = X , with 77 given by 
(10), dQ = J7 _1 dP and X given by (5). We note that 77 G D°°> 2 (see, e.g., [9] Lemma 
6.3.1 and [2] for the expression of Malliavin derivatives via the transfer principle). We 
can particularly deduce that each random variable Vfc , recursively defined by Vq = T) and 
V k+1 = (DV k , l[ ,t])si ^ k > 0, belongs to B 1,2 . 

In (22), the deterministic terms f* +h \v - u\ 2H ~ 2 dv, (t + h) 2H - t 2H and 

t+h 

((v - t) 2 "- 1 - v 211 " 1 )^ 211 + v 2H -(v- t) 2H ) dv 

have a Taylor expansion in h of the type (15). Lemma 4 allows the first term of the 
asymptotic expansion to be obtained using the fact that J^ +h (j>( s ) ds = htj)(t) + o(h) for 
any continuous function <j>. By a recursive argument, again using Lemma 4, we finally 
deduce that (15) holds. This concludes the proof of Theorem 2. □ 
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